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Peculiarity of the Kondo effect in CeajLai-^Be is investi- 
gated on the basis of the scaling equations up to third order. 
For the case where the J 1 -/ 2 charge fluctuation enters in ad- 
dition to the f 1 -/ one, the effective exchange interaction 
becomes anisotropic with respect to the orbital pseudospins 
which represent the two different orbitals in the T$ ground 
state. Because of different characteristic energies for electric 
and magnetic tensors, scaling with the single Kondo temper- 
ature does not apply to physical quantities such as the re- 
sistivity and magnetic susceptibility. Possibility of a bizzare 
phase is pointed out where the RKKY interaction leads to the 
spin ordering without orbital ordering. This phase serves as 
a candidate of the phase IV which is observed to be isotropic 
magnetically. 



I. INTRODUCTION 

Orbital dynamics in a number of heavy-fermion sys- 
tems has recently attracted great interest jj] || . In the 
presence of orbital degeneracy, orbital (electric) tensors 
in addition to spin ones have an opportunity to be active, 
and shows rich phenomena resulting from entanglement 
with spin (magnetic) degrees of freedom. The orbital 
degrees of freedom couple with lattice ones, and often 
lead to the Jahn- Teller effect. In other quench- 
ing mechanism for these degrees of freedom such as the 
orbital Kondo effect prevents static distortion of the lat- 
tice. Hence study of the coupled spin and orbital fluc- 
tuations should be important for understanding heavy- 
fermion systems as a whole. 

Under a high symmetry such as the cubic one, there 
are chances for the orbital degeneracy to remain. A typi- 
cal example is the cubic compound CeBg : its crystal- held 
ground state is the Tg quartet which consists of two de- 
generate Kramers doublets. The excited doublet T 7 is 
well separated by about 540K and plays little role in 
low-energy physics. CeBg exhibits curious phase diagram 
at low temperature in magnetic held | ^— [Tc|] . The phase 
boundary between the paramagnetic phase, called I, and 
the antiferro-quadrupolar phase called II shows unusual 
dependence on magnetic held: transition temperature in- 
creases as magnetic held increases. The magnetic held 
dependence of the phase II has been ascribed either to 
the intersite interactions between higher-order multipolcs 
Jil"|,fl| , or the quadrupolar fluctuations ||^| . 

In the course of systematic dilution study jij] |l4| , 
a strange phase, called IV, is observed recently in 



Ce^Lai-^Bg. The magnetic susceptibility shows a cusp 
on entering the phase IV from the paramagnetic phase 
I with decreasing temperature. This suggests that the 
Neel state is present here. In contrast to the phase III 
which has both antiferromagnetic and quadrupolar or- 
ders, the phase IV has a very small magnetic anisotropy 
in the susceptibility |lg| ] and almost no magneto resis- 
tance [ pd[ . In realizing this phase the interplay between 
intersite correlation and on-site Kondo effect seems to be 
essential. 

In a previous paper jDsf , one of the present authors 
noted the importance of the orbital Kondo effect in un- 
derstanding the phase IV; the orbital Kondo effect can be 
active even in the presence of a spin ordering. Then a new 
method was proposed to perform perturbative renormal- 
ization to arbitrary higher orders. In order to quantify 
the idea, however, one must estimate actual magnitude of 
relevant interactions under the realistic point-group sym- 
metry. Up to the present, investigation of the Kondo- 
type interaction in the presence of crystalline-electric- 
ficld (CEF) effects has been made by various methods 
|]l|,p| jlqjl7|] . It has turned out essential to take into ac- 
count the splitting of the localized states in accordance 
with the point-group symmetry As a result, the ef- 
fective exchange interaction needs many parameters for 
characterization. In the conventional held theoretical 
scaling procedure, it is tedious to deal with such large 
number of parameters. 

In this paper we generalize the scaling method of Ref. 
Jl5| so as to be applicable to arbitrary point-group sym- 
metry. The scaling equations for symmetry adapted cou- 
pling constants are written down generally in terms of 
structure constants of the relevant Lie algebra. As a spe- 
cific case, the third-order scaling for the cubic symmetry 
is performed explicitly with minimum amount of interme- 
diate steps. We show that the pscudospin representing 
the orbital moment has an anisotropic exchange interac- 
tion as a result of scaling. On the other hand the mag- 
netic pseudospin remains isotropic in the case where the 
J 1 -/ fluctuation dominates the Z 1 -/ 2 one. Since the 
third-order scaling does not work for a strong-coupling 
fixed point, the local Fermi liquid cannot be identified 
by the present approach as it stands. However, with 
available knowledge from various sources, we can almost 
certainly classify all the fixed points of the model. It 
turns out that there is an unexpected symmetry in the 
exchange interaction which interchanges both the spin 
and orbital indices simultaneously. We show that this 
hidden symmetry is specific to the SU(2)xSU(2) sym- 
metry which is relevant to the Tg CEF state. 
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This paper is organized as follows: In the following 
section, we describe the method to derive the effective 
exchange interaction starting form the Anderson model. 
This section is mainly adaptation of earlier treatment fl(f | 
to respect the point group symmetry from the beginning. 
The details of explicit derivation are given in Appendix 
A. In section III we derive the third-order scaling equa- 
tion in the most general case of the point group. In sec- 
tion IV we apply the general result to the specific case 
of the Ce x Lai_ x B6 system and give detailed analysis of 
the scaling equations. The final section summarizes the 
present paper with implication to identifying the myste- 
rious phase IV. 



II. EXCHANGE INTERACTION IN 
IRREDUCIBLE REPRESENTATION 

The atomic structure of a single magnetic ion with / 
electrons is treated by the CEF theory together with the 
Russell-Saunders {LS) coupling scheme. In this scheme 
the local state is specified uniquely by the following quan- 
tum numbers: the number of / electrons /", the orbital, 
the spin and the total angular moments L, S, J, and the 
irreducible representation T (abbreviated as irrep) of the 
double point-group with time-reversal operation and its 
component 7 together with the branching multiplicity 
label Q, which is required for the case where the ir- 
rep r occurs more than once (e.g. -D5/2 = I^ ffi 2I\ 
for tetragonal symmetry). We often abbreviate the lo- 
cal states \ f n LSJQT^f) as \4>nl) for notational simplic- 
ity. In the same manner the one-particle state of / elec- 
tron is described as |£A) where £ is the abbreviation of 
(£ = 3, s = 1/2, jqA) with q being the multiplicity label, 
and A specifying the irrep of the one-particle state. The 
corresponding creation operator is written as /^ A . It is 
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useful to express a Bloch state conduction operator c 
using c k £ X with the symmetry adapted basis around the 
impurity 

The Anderson Hamiltonian then takes the form: 



H = Hk+ Hf + -ffhyb, 



k£\ 



0n7 



#hyb = [ V K C U\f^ + h -' 



(1) 

(2) 
(3) 
(4) 



where we restrict the hybridization between the local and 
the conduction electrons to the same set of symmetry 
(not only A c — Af but also £ c = £/). We regard the 
label £ as the channel index of independent scattering 
processes and A as the internal degrees of freedom which 
is responsible for the Kondo effect. 



If the most stable configuration /" is well separated 
from f n±1 ones, one can restrict the model space E]| 
to the multiplet (f n LSJ) by integrating out the virtual 
charge fluctuations to the f n±1 configurations. The ef- 
fective Hamiltonian in the model space is written as 



H cB = H k + VH f V + H cx + G(y 4 ), 
H cx = \a)(a\V(E„- Hf^QVlb)^, 



(5) 
(6) 



a,b£M 



where V denotes the projection operator to the model 
space M and Q = 1 — V . We have neglected H k in 
the resolvent of intermediate state assuming that the 
relevant conduction-electron states have energies smaller 
than those for excited states of if/. The exchange pro- 
cess via the excited configurations is shown schematically 
in Fig.|l|. Then the exchange interaction generally takes 
the form 



J ia cVc- |0')X Q < 



(7) 



where meaning of the indices i and a are to be specified 
later. The matrices x* and X" describe the transition pro- 
cesses for conduction electrons and the localized states, 
respectively. Summation over quantum numbers of each 
state is implied by the matrix multiplication. The num- 
bers of independent matrices are 14 2 — 1 for the conduc- 
tion electron and d j — 1 for the multiplet (f n LSJ) with 
the degeneracy dj. The unit matrices which give the po- 
tential scattering do not enter the exchange Hamiltonian. 

The invariance under the point-group operation is im- 
plicit in eq. (0) . To make it explicit Jl6[ | , we introduce the 
irreducible tensor operators M] defined as follows: 



k'X'kX V 



X { ll(<t>'n<f>n)=YWnl'){<f>nl 



7'7 



A'* 
A'* 

■p/* 
/* 
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A* 
6* 



A r 

5 7 



(8) 



where the one-column and the three-column brackets in- 
dicate 2jm and 3jm symbols with the multiplicity r or 
t jfJJ. These are natural extensions of the lj and 3j 
symbols of Wigner. It is clear that the tensor operator 



(r) 
V A*S 



(£'£) transforms like the ket \rA*S*), or equiva- 
lently like the bra (rAS\, under the point-group opera- 
tions. The irreps of the localized-state tensors are de- 
termined so as to have the finite 3jm symbols. Namely, 
decompositions of the direct product T'* ® T contains A. 
Similar decomposition also determines the irreps of the 
conduction-electron tensors. We write the exchange in- 
teraction Jia as another matrix g in the invariant form. 
Namely cq.(0) is equivalently written as 

£' ' (,4>' n <t>n Art 

' A 



><E S (9) 
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Here the 2jm symbol is inserted in the exchange inter- 
action. With the rotational symmetry, eq.(||) reduces to 
the known result involving the spherical tensors 16 3]. 

The magnitude of the coupling constants, which will 
be obtained explicitly in Appendix A, are given by 

= £ Ai rt) (0 n+1 )(^n/ e ii^ +1 )<0„ +1 |i4ii0„) 

+ £ ^(^-iX^II/Jll^-iX^-ill/e'HW, (io) 



where A^\(j) n ±i) is of the order of \V\ 2 /[Ef(<j) n ±i) — 
Ef((j) n )]. The first summation is taken for the excited 
configuration f n+1 such that both the direct products 
r' Cg) A' ® r + and T + ® A F contain the identity repre- 
sentation. Similar selection rule is available to the 
configuration po|j . 

Hereafter we restrict our discussion for simplicity to 
the case where the irrep AS is real and the correspond- 
ing 2jm symbol can be set to unity. Moreover, both 
the initial and the final states belong to the ground-state 
CEF multiplet 4> g . Then the localized-state tensors be- 
come Hermitian and all have definite signs with respect 
to the time-reversal operation. The conduction-electron 
tensor does not have the definite sign in general, since the 
time-reversal operation interchanges £' with £. However, 
if one considers the matrix of the tensors in the combined 
space £' £ which are defined as 



(£' = 



*S(£'0 + 41(^0 (otherwise) 



(11) 

then the matrix of the tensors have definite signs under 
the time reversal [fl6f . In the above discussion we use 
the fact that the coupling constants are real and sym- 
metric against interchange of £' and £ which is ensured 
by the hermiticity and the time-reversal symmetry of the 
exchange interaction. In the treatment above the same 
coupling constant is automatically imposed for the dif- 
ferent channels. This property was essential to derive 
the two-channel Kondo Hamiltonian for tetragonal and 
hexagonal symmetries pcfl . We rewrite the exchange in- 
teraction in the restricted case as 



41(^)4Sou (12) 



(£'£) Art 



where the summation of channels is taken for combina- 
tion of a pair of £' and £ because of the symmetrized 
expression introduced by eq.(0). 



III. SCALING EQUATIONS FOR GENERALIZED 
EXCHANGE INTERACTION 

We derive the scaling equations up to third order for 
the exchange interaction in terms of the irreducible ten- 
sors. The index i in eq.(^) is abbreviation of (rAS, 
and the index a of (tA5, <fr g ). In dealing with the ma- 
trices representing the irreducible tensor operators, we 
note that the matrix x* actually depends only on a sub- 
set of quantum numbers (A, A', A), and the matrix X Q 
on (A,r s ) where T g is the irrep of <j> g . We write the di- 
mension d(A, A', A) of x l simply as di and the dimension 
D(A,r g ) of X Q as D a . To treat all matrices on an equal 
footing, each matrix should be embedded in a space of 
14 x 14 matrix. The exchange Ji a in eq.(0) is defined as 

Jia = Sa l° r (Ac<5 c ) = (AfSf) and zero other- 

wise. Here Ji a does not depend on the component index 
S. 

The matrices satisfy the orthogonality relation: 



Tr(xV) =di5ij, 
Tr(X a X^) = D a 5 a0 . 
The commutation rule is given by 

[ xI > xJ ] = *E j-/«*x*> 

k k 



(13) 

(14) 
(15) 



Equivalently the structure constant is given explicitly by 

/^--zTr([x\xV), (16) 
F aM = -iTr([X a ,X^). (17) 

It is obvious that the structure constants fijk and F a py 
are completely antisymmetric against interchange of a 
pair of indices. 

According to the renormalization formalism based on 
the open-shell Rayleigh-Schrodinger perturbation theory 
Jl5]| , a change of the band cut-off E c induces the following 
expansion of the effective interaction matrix: 



(2) 



Sh 



(3) 



'i" 1 u "int "f u "int ' " ' (1^) 



where the superscript indicates the order of the bare cou- 
pling constant, and the lowest-order matrix h ox is given 
by 



hex — ^ ^ Jiot^ X . 



(19) 



The second-order contribution shown in Figj^(a) and 
^(b) to the effective interaction is given by 

/<5(ln E c ) = E J *° J ^ ^ > xJ l ( XQX/3 ) 

ij a/3 

^^EE J «^[x l ,x^][X«,X' 3 ], (20) 

ij a(3 
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where we have interchanged the dummy indices (i, a) «-> 
(j, (3) in deriving the second equality. 

Figure ||(c) shows a diagram for the third-order con- 
tribution to the effective interaction. One should also 
take into account the "folded diagram" |l9] shown in 
Fig.^(d). Taking average of the original and Hermite- 
conjugate counterpart of the folded diagrams, we obtain 

6h^/5(hxE c ) = \ J2 E J - J ^ J fc7 x l Tr(x^x fe ) 

ijk af3y 

x {X^X^X"] - [X^X^X 7 } 
= \ E E ^^-/3 Jfc 7 x l Tr(x^x fc )[X' 3 , [XT, X"]]. (21 

ijk uf3f 



Here we have also interchanged the dummy indices 
«-> (fc,7) in deriving the second equality. 
By computing the commutators we find that the hj nt 
has the same matrix structure as the lowest order part 
h cx . Thus, we obtain the scaling equations in terms of 
the structure constants as 



d 



J - B {2) + /3 (3) 

■J at — Mae ~ i J at I 



ij a/3 



1 



Pa} ~ r. r-. Jga ^""^ djJjfiJj^ ^ j n _F al0 F e ps, 



2D e ^ 



s Ds 



(22) 



with £ = In 2? c . This is the most general form of the third- 
order scaling which is valid for any point-group symme- 
try. We emphasize that non-commuting property of ten- 
sor operators is concisely taken into account in terms of 
structure constants of the underlying Lie algebra. 



IV. APPLICATION TO CUBIC SYMMETRY 

In this section the scaling analysis developed in the 
previous sections is applied to the case of dilute system 
Ce K Lai_a;B6 which exhibits the remarkable entanglement 
of magnetic and electric tensors in static and dynamic 
properties. We first derive the exchange interaction with 
use of pseudospins and then discuss nature of scaling. 



A. Exchange Interaction 

The magnetic ion Ce 3+ (J 1 ) lies in the cubic-symmetry 
(Oh) CEF. The degeneracy of the ground multiplct 2 F 5 / 2 
is lifted to the excited doublet and the ground quartet 
r 8 . The partial waves of conduction electron are also 
classified by the cubic symmetry: 



F>5/2 — r? 
D7/2 = F 6 



r 8 , 



r«. 



(23) 
(24) 



We define the basis sets for both conduction (j) and lo- 
calized (J) states as follows: 

(i) j, J = 5/2 



|r 7 : (T,i)> = V 6 
|r 8 : (+!,+ !)> = 
|r 8 : (-!,-!)) = 



± 5 - 

2 
5 

1 

±- 



T 2 



±- 



) (ii) j = 7/2 



|r 6 :(U)> = ± 
|r 7 :(U)> = ± 




|r 8 : (+!, + !)) =± 



(25) 
(26) 
(27) 

(28) 
(29) 
(30) 
(31) 



where irrelevant quantum numbers have been omitted. 
To specify the components of irrep, we have used the 
symbols j and j for the time-reversal partner and the ex- 
tra orbital labels ± for the Tg irrep. The relative phases 
of the basis are chosen so that the Kramers pair trans- 
forms like the spin 1/2 under the time-reversal operation 
9, Le.,9\ T) = |i>,0| I> = -It>- 

Before we express the irreducible tensors, it is conve- 
nient to introduce two pseudospins a a and r l (each of 
them is defined as usual convention of Pauli matrices) 
which act on the Kramers and the non-Kramers pairs, 
respectively, without changing the other degrees of free- 
dom The six pseudospin operators are classi- 
fied by the time-reversal operation. Due to the definition 
of the pseudospin and the basis, it is easy to obtain their 
property under the time reversal: 



odd (magnetic) 
even (electric) : 



: a 



a v ,a z ,T y , 



(32) 



Notice that the pure imaginary operator r v has the same 
transformation property as magnetic moment. 

Once we fix the basis sets, we can express the irre- 
ducible tensors (lj]]2ll] in a concise way by using the pseu- 
dospins. We define cP and r° as unit matrices in the spin 
and the orbital spaces, respectively, and the linear com- 
binations of t x and t z as 



(33) 
(34) 



All tensors for possible combination of the basis sets are 
given as follows: 
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(i) r 6 x r 6 

r 4m : [a x ,a y ,a z ], 

(ii) r 7 x r 7 

r 4m : [a x ,a y ,a z ], 

(iii) r 8 x r 8 



T2m 

r 3e 



rV° ] , 



4 m 
(2) 



A 4m 
I^m 
I^e 



(iv) r 6 x r 7 



2 ( 



5 m. 






a 


a u 










a x 











O y 






(35) 



(36) 



(37) 






cr 2 


a z 






(38) 






a y 























( o 


a z 





- 












V° 









5m 






-V3V 


a x 















, 






\fia y 


a y 










a y 









' 





a 2 











> az 









. (40) 



where the subscript m represents magnetic and e does 
electric. 

The exchange interaction can be cast into the form 
cq. (^) by using the explicit form of the irreducible ten- 
sors. The number of independent coupling constants 
for each localized-state tensors is evaluated by counting 
possible combination of scattering channels in the same 
irrep. Precisely, there are 56 in total: n^m — ?^5e = 3, 

(2) 



n- )e = 8, n 5m = 10 and n^ n = ri^ = 16 



We note that f 1 -/ charge fluctuation involves only 
the scattering channels of (j — 5/2, Tg) symmetry in 
the exchange interaction because of the selection rule, 
eq.([l(]). Thus, if one ignores f 1 -/ 2 charge fluctuation, 
the exchange interaction reduces to that of the SU(4) 
Coqblin-Schrieffer(CS) model. 



(v) r 6 x r 8 
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3<? 



4m 



5m 
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V3a y 


a y 










<T y 












a y 


-Vs<7 y 


a y 








-V3a y 









(vi) r 7 x r 8 






-V3a x 


a x 


-^JZa x 








a x 





t 





i o 





a z 



























a x 


V3a x s 


a x 








VZa x 





o , 





/» 


a z 

















Vo 









(39) 



3e 
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a x 


V3a x " 


a x 








V3<j x 
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B. Explicit Form of Scaling Equations 

In deriving the scaling equations explicitly, we restrict 
ourselves for simplicity to the case where the f 1 -/ charge 
fluctuation dominates over the Z 1 -/ 2 one. In this case, 
the Z 1 -/ 2 fluctuation gives two different corrections to 
the SU(4) CS model: (i) modification of coupling con- 
stants breaking the SU(4) symmetry, and (ii) generation 
of additional scattering channels besides (j = 5/2, Tg). 
The latter correction gives very low characteristic energy 
as compared with the one given by the J 1 -/ fluctua- 
tion and hardly affect the renormalization of the SU(4) 
CS model. On the other hand, the former correction 
changes the renormalization flow qualitatively and gives 
rise to multiple characteristic energies even though it is 
small. 

We take into account the effect of J 1 -/ 2 fluctuations 
only for such processes that are absent in the SU(4) CS 
model. Then there appear the exchange interaction with 
the irreps A = 2m, 3e, 5m, 5e, 4m by the decomposition 
of Tg (g> Tg. Of these the last one A = 4m has the mul- 

(rt) 

tiplicity 2 which we distinguish by using the matrix g\ 



with r, i = 1, 
interaction: 



2. We deal with the following exchange 



H cx — — 



9a(xX) a + 53 £ \ x rX t ) 

A=£4m r,t=l 



4 m 



, (41) 



where we define the summation of the components as 



5 



(x r X t )i 



s 



(42) 



with the scattering channel, £' = £; = (j = 5/2, Tg) and 
the local ground-state configuration <f> g = (J 1 , 2 F 5 / 2 , Tg). 
It is noted that the SU(4) CS model is reproduced by 
setting all coupling constants as equal to g except for 

(12) (21) ,_ ^ 

We use the formulas (|16|). ( J17| ) and (22) together with 
the definition of tensors (|37|) to obtain the set of scaling 
equations: 







dl 

d 1 

^53e — -^^ffamffS 

9_ (ii) 

»54m 



l 2 , o (22) n 1 ro 2 . Q / ( 12 )2 . (22)2 2 s. 



■ 355e(54m ) + 55m)] + ^93e[2(gl m + 53e) + 3(54™^ 



9ir 



(22)2 



2 

9bn 



2<&)], 



dt 

d_ (12) 

<9 £ 54m 



d_ (21) 

m 9im 
9_ (22) 

d£ 9im 



m 95r, 



[45im )2 



(22)2 
54m 



2(5L 2)2 



654m ^95r, 



r (11)3 _ (12) (21) (22) (11), (21)2 
li/4m 54m f4m 54m T 54m V54m 



1 2 
i + 55m 

(22)2 
54m 



b 4 5 | e ] 
5eL 



55m)]: 



"2^"4m U/4m — <Jy5mJ — z 54m 54m J 

4. Ir„( 12 ) 3 _4„( n )„( 21 )„(22) 

+ ^[9 Am 4 54m 54m 54m 

l \n (12 ^(J 2 ^ In -\ 9«( 1:L ^( 21 )l , [„(21)3 (11) (12) (22) 

'77l54m 154m ~ d 55mJ ~ ^54m 54m J [54m — 54m 54m 54m 



(12) r 4 (21)2 
4m I ^J4m 



+ 54 



(22)2 



2(52 2 m 



5 3 2 e 



55e)+ 55 2 m}] J 



(43) 
(44) 

(45) 

(46) 



9 4 



(21)/„(H)2 
\94m 



5 4 



L 2)2 +5 5 2 e + 5 5 2 J], (47) 



,(12)„(21) 



2 [54m' 5 I,, 



„(H)/„(22) 
54m \54m 



J[5i 2 m 2) {45im )2 



(22)3 



±rVo (11) o (22) 

oL*H»4m 54m 



(12) (21)s 
54m 54m ) 



T[55m{4(54m )2 



^55e = -[2g { l^g be + g 3e (g 4m 



(22) 



' 55m)] 



355m) + 2(53e55e + 52m55m)] 

54m + 2(52m + 5 3 e + 5 5 e) + 5 5 5 m} + 5 4 r, 

of (22) 1 \ , (11) 1 

" Z \9im 92m + 53e55eJ + 54m 55mJ 

/ (12)2 . (22)2x . „, 2 1 2 , 2 \i 
(54m +54m ) + 2 {92m + 93e + 95e)\ 

■i , „ JofnC 11 ) 2 J. „( 21 ) 2 , „2 \ , „(12)2 (22)2 2 1, 
55e|2(54m + 54m +53e)+54m + 54m +55m/J' 



(21)2 
54m 



^54m 54m 54m J' 



' 55m] J 



(48) 

(49) 
(50) 



We note that these apparently complicated expressions 
follow straightforwardly from eq.(22). The correctness of 
the expression has been checked by taking various limit- 
ing cases. 



C. Nature of scaling 

Let us discuss implication of the scaling equations. 
The set of equations has eight stable as well as saddle- 
point fixed points, which are summarized in Table Q. Be- 
cause all fixed points have the relations g^* — g% m and 

9im* — 54m = 0' the exchange interaction at fixed 
points reads 



where the coupling constants have been redefined as fol- 
lows: 



54m ) = J > 93e = Kj_, g 2m = K z , 



J22) 
54m 



55r, 



55 



I,. 



(52) 



1 



4 ^ 



4 

T X T X 



c f y 



<T?<T?, (51) 



The pseudospins for conduction electrons are rewritten 
explicitly as 

± U 

T c a c = C fc' m' a' Pm' mPa' a c km<7 , 

k'k m'm a' a 

(i,a = 0,x,y,z), (53) 

where p a with a = x, y or z denote the Pauli matrices 
and p° does the unit matrix. Similar definition is also 
used for the pseudospins for localized states. In the rede- 
fined expression of exchange interaction (|5l|), difference 
in time-reversal characters of the irreducible tensors ap- 
pear as the exchange anisotropy of pseudospins. This 
kind of cr-r double tensor exchange model has been stud- 
ied in the literature. However unnecessary imposition on 



G 



the parameters caused ambiguous conclusion about the 
fixed points fj2f . 

The two fixed points in a given row in Table || are es- 
sentially the same with each other since the one with the 
upper sign changes into another with the lower sign by 
a unitary transformation which changes simultaneously 
the signs of the transverse couplings, K± and Jj_. The 
groups with and without a prime are related to each other 
by the transformation I z ) <-> — (I±, I z ). This hidden 
symmetry will be discussed later. 

In the absence of the last term in eq. ( ]5l| ) , the a and r 
spaces are decoupled. Thus the fixed points (i), (ii) and 
(iii) correspond to: (i) the non-Fcrmi-liquid (NFL) fixed 
point in a space, and (ii) the Ising and (iii) the NFL fixed 
points in r space, respectively. 

The last term couples a and r spaces. In the absence of 
K± , the term of I z leads to the fixed point (vi) where the 
NFL appears only in a space with the coupling constants 
Jc« = J + Iz T c T f ! while the term of I± to the fixed point 
(vii) . In the presence of I z or I± together with the trans- 
verse coupling K±, the coupling constants flow to the 
fixed point (viii) with the SU(4) symmetry. It is known 
that the Coqblin-Schrieffer model does not have the NFL 
fixed point. Then, the finite magnitude of the fixed-point 
coupling is an artifact of the third-order scaling. The cor- 
rect fixed point should be at J* = K± = K* = oo and 
gives the local Fermi liquid. The stability of the saddle- 
point fixed points against each type of perturbation is 
summarized in Table [n]. 

Let us discuss renormalization evolution and charac- 
teristic energies for some simplified cases. First, we con- 
sider the isotropic case: K± = K z = K and 7j_ = I z = I. 
In this case, the scaling equations are reduced to those 
discussed in Rcf. [|15| as 



d_ 

di 
d_ 

di 
d_ 

W 



j 



K 



-a-j)(j 2 +3i 2 ), 

-(1- K)(K 2 + 3I 2 ), 

-21 (J + K) + I(J 2 + K 2 + 2I 2 ) 



(54) 
(55) 
(56) 



In the absence of the coupling I, it is obvious that there 
exist two characteristic energies corresponding to each 
Kondo effects: T^/E c = Jexp(-1/J) and Tg/E c = 
K exp(—l/K). On the other hand, the case of SU(4) 
symmetry with I = J = K gives single characteristic en- 
ergy, Tg s /E c = I 1 / 4 exp(-l/47). It is noted that T£ s 
is larger than T K ' since the number of the screening 
channels is larger. The renormalization evolution of the 
coupling constants is shown in FigJ| for the initial cou- 
pling constants Jo = 0.2 and Kq = 0.1 with Iq = 10~ 6 , 
10~ 4 and 10 -2 . The three different characteristic ener- 
gies corresponding to J, K and I merge together as Iq 
increases. The average of the characteristic energies also 
increases as Iq increases. 

Next, we focus on the difference in time-reversal prop- 
erty and put J m = J = K z = I± and J e = K± = I z . 



The scaling equations are simplified in this case as 

§i J ™ = -( 1 - Jm)(3Jl + J 2 e ), (57) 

^J e = -2J e [2J m - {J 2 m + J 2 e )]. (58) 

In the absence of the electric coupling J e , the characteris- 
tic energy corresponding to the magnetic Kondo effect is 
given by T^/E c — Jm 3 exp(— l/3J m ). T^ n is also smaller 
than T^? s since only the magnetic channels contribute to 
the Kondo effect. The renormalization evolution is shown 
in Fig.[|for the initial coupling constants J m o = 0.05 with 
various values of the ratio a = J e o / J m o ■ The two differ- 
ent characteristic energies corresponding to J m and J e 
merge together as a approaches unity. 

Before closing this section, we discuss the hidden sym- 
metry / «-> —I. Let us consider the following particle- 
hole transformation for the conduction electrons: 



Ckr, 



-kia ' 



(59) 



where the matrices S and T are unitary. The transfor- 
mation does not change the kinetic-energy part of the 
Hamiltonian unless the conduction band is asymmetric 
with respect to the Fermi level. On the other hand, the 
spin operators er* and t" are transformed respectively as 



T c ~ ¥ c k'm't?Pm'm c kmcn 
k'k am'm 

a _ _ \ " \ " f -q 
°e "~ * 2^1 2-^1 C k'ma'Pcr'cr c kma, 
k'k ma' a 

T c CJ c > ~ ^ ] y ] c k'm'cr'Prn'mPa'a c km<j, 



(60) 
(61) 
(62) 



k'k 



m ma a 



where the minus signs come from the anticommutation 
of conduction-electron operators. The transformed spin 
operators with bars are given by 



(63) 



Let us assume the presence of unitary transformations S 
and T which make the transformed spin operators pro- 
portional to the original one. Namely, we have p l = Ap l 
and p a — Bp a for all components where A and B are 
independent of the component indices. Then the trans- 
formed exchange interaction is equivalent to the original 
one with the coupling constants J = — AJ, K = —BK 
and I = —ABI. There indeed exist such transforma- 
tions in the case of the SU(2) symmetry; the matrices 
are given by S = T = p y which give A = B = — 1. 
However, in the case of SU(7V) symmetry with N > 2, 
such transformations cannot be found since the rank of 
the SU(.ZV) symmetry is higher than 2. We note that the 
hidden symmetry for SU(2) plays a special role to reject 
the NFL fixed point which is present in the large limits 
of both N and M in the SU(iV)xSU(M) exchange model 
discussed in Ref. fn5|. 
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V. SUMMARY 

In this paper we have proposed a systematic proce- 
dure to derive the exchange interaction and the scaling 
equations up to third order with maximum use of the 
point-group symmetry. Various moments are described 
by the irreducible tensors with proper time-reversal prop- 
erty. Then the exchange interaction is obtained in the 
scalar- product form of the irreducible tensors. The third- 
order scaling equations for symmetry adapted coupling 
constants are written down generally in terms of struc- 
ture constants of the relevalnt Lie algebra. 

The procedure is applied to the case of Ce^Lai-^Bg 
where the CEF ground quartet Tg shows rich phenomena 
resulting from the entanglement of electric and magnetic 
tensors. We derive the exchange interaction integrating 
out the charge fluctuations to /° or f 2 configurations. 
We discuss the nature of scaling for the case where the 
J 1 -/ fluctuation dominates over the J 1 -/ 2 one. In this 
case, main contribution to the Kondo screening comes 
from the partial waves with j = 5/2, Tg symmetry of 
conduction electrons. 

As a result of scaling, the effective exchange interaction 
is described by the cr-r double exchange model with the 
exchange anisotropy in r space. The anisotropy comes 
from the fact that one of three components of the orbital 
pseudospin t v has different time-reversal property from 
others. The effective exchange model has non-Fermi- 
liquid fixed points in the absence of the cr-r coupled term. 
In the presence of the coupled term the system flows to 
the local Fermi-liquid fixed point. The different charac- 
teristic energies corresponding to the exchange couplings 
for either a or and r rapidly merge together with increas- 
ing magnitude of the coupling. If one considers only the 
time-reversal difference, the system has also two different 
characteristic energies corresponding to the electric and 
the magnetic exchange interactions. We identify the spe- 
cial symmetry of SU(2), which plays an important role 
to reject the non-Fermi-liquid fixed point resulting from 
the argument of SU(iV)xSU(M ) exchange model in the 
large limit of N and M |l5| . 

Concerning the physical quantities in the intermedi- 
ate regime toward the fixed point, magnetic and electric 
quantities can have different energy scales. Therefore, 
conventional argument with the single Kondo temper- 
ature has to be revised in the presence of the orbital 
degeneracy. This should affect previous analysis of the 
resistivity and magnetic susceptibility. For this purpose, 
quantitative argument is required beyond the perturba- 
tive scaling analysis. In the forthcoming paper we use the 
numerical renormalization group |23| ] to derive physical 
quantities at finite temperature. 

In the presence of intersite interactions, the orbital 
and spin degeneracies can be split by a molecular field 
associated with the long-range order. It is known that 
Ce^Lai-^Bg has at least four different phases. Of 
these the phase II is characterized by the orbital (or 



quadrupole) order without magnetic order. The mag- 
netic order present in the phase III seems to accompany 
the orbital order. If the molecular field corresponds to 
the internal magnetic field, not only spin but orbital de- 
generacy are lifted because Zeeman splitting depends on 
each orbital. Provided that the intersite spin exchange 
interaction is stronger than the orbital one, we can think 
of the situation where the spin Kondo effect gives way 
to the long-range order but the orbital Kondo effect does 
not. In such a situation the magnetically ordered phase 
should keep the nearly isotropic magnetic property as 
in the paramagnetic phase. We expect that a detailed 
quantitative study with the scenario described above will 
provide us with understanding the nature of the curious 
phase IV in Ce^Lai_ x B6. 
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APPENDIX A: EXPLICIT DERIVATION OF 
EXCHANGE INTERACTION 

Here we derive the explicit expression of the coupling 
constants by applying the second-order perturbation to 
the Anderson model (|Ij). The exchange interaction via 
the n + I configuration is given by 

H ^ = - E E 7 + E 4rv c ^i^y)(0n7i 

with 



E f (4> n+1 ) - E f {<j> n ) 

x <0:ji/eii^ + i)(0„ + iii4ii^)- (A2) 

The reduced matrix elements of the creation operator 
can be calculated by using e.g. the Racah factorization 
lemma in addition to the coefficient of fractional parent- 
age (T9|2|. 

The summation of the products of the 3jm symbols in 
eq.( |Al[ ) can be converted to the summation of another 
combinations of 3jm symbols which relate to the prod- 
uct of the irreducible tensors |l6| . For this purpose the 
following identity is useful: 
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E 

7+ 

E 



T+ A' i 
A* A'* A 



A r 

7* A 7 



r+ r r' , 

ASrt *• *? 9 r 



|A|{riArg}{r*A*r't}{r*} x 



A' 
A' 



A'* A* A 
A'* 5* A 



r'* a r 

(A3) 



where the 6j symbol in the point-group irrep has ap- 
peared together with the 2j and 3j phases 19 . By using 
the definitions of the irreducible tensors, eq.(j|), we ob- 
tain the coupling constants for the n + 1 process as 



sit = E'+{ 



A* A'* A 

r+ r r' 



|A|(-1) P , (A4) 



tq' qr 



<Pn + l 

with the phase 

(-i) p = -{riAr g }{r*A*r'i}{r*}. 



(A5) 



Similar argument for the n — 1 configuration provides 
the coupling constants for the n — 1 process: 



^ = E 7 - 



A* A A'* 

r_ r r 



|A|(-1) 



F" 



(A6) 



tqq'r 



with 



x (^ii4n^-i)(^-iii/^ii^), 



(-i) p ' = {A*}{A*}{r*} 



{r_r'*A 9 }{A*AAV}{r'*Art}. 



(A8) 



Thus, the full form of the exchange interaction is given 
by 

H ™ = E Esa^'^") 
' A 



E 7 (A9) 
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FIG. 1. The exchange process via the excited configura- 
tions <^>n±i7±- The solid (dashed) line denotes the conduction 
(one-particle /) electron. The double solid line represents the 
localized state in the projected space. 



(a) (b) 




FIG. 2. Scattering processes (a), (b) in second order and 
(c), (d) in third order. The solid line shows a conduc- 
tion-electron state, while the dashed line the local electron. 
The backward propagation of the dashed line in (d) is char- 
acteristic of the folded diagram. Each index denotes the 
corresponding coupling constant. The intermediate conduc- 
tion-electron states is required to have energies near the 
cut-off E c . 




FIG. 4. The renormalization evolution of coupling con- 
stants for the initial magnetic coupling constant J m o = 0.05 
with increasing the ratio a = J e o/JmO, Je being the electric 
coupling constant. The two different characteristic energies 
corresponding to J m and J e merge together as a approaches 
unity. 

TABLE I. The stable and the saddle-point fixed points for 
scaling equations. The coupling constants in the redefined 
notation are written in the second line. 

(11)* * « (22). _ * * (12)* _ (21)* 

54m 9Ze 92m 9 im — 95m <?5e 9im ~ 9 4m 





J 


K± 


K z 


I± 


/* 




(i) 


1 

















(ii) 








K z0 











(iii) 





±1 


1 











(iv) 


1 





K z0 











(v) 


1 


±1 


1 











(vi) 


1 





K z0 





+1 





(vi') 


I 





K z0 





-1 





(vii) 


I 





1 


±1 








(viii) 


1 


±1 


1 


±1 


+1 





(viii') 


1 


±1 


1 




-1 







TABLE II. The destinations of the saddle-point fixed 
points (i)-(vi') in the first column against each type of per- 
turbation. 





J > 


J < 


K x 


K z 


Ix 


I z > 


h < 


(i) 






(v) 


(iv) 


(vii) 


(vi) 


(vi') 


(ii) 


(iv) 


(") 


(iii) 




(vii) 


(vi) 


(vi') 


(iii) 


(v) 


(iii) 






(viii) 


(viii) 


(viii') 


(iv) 






(v) 




(vii) 


(vi) 


(vi') 


(v) 










(viii) 


(viii) 


(viii') 


(vi) 






(viii) 




(viii) 






(vi') 






(viii') 




(viii') 






(vii) 






(viii) 






(viii) 


(viii') 



FIG. 3. The renormalization evolution of coupling con- 
stants for the initial coupling constants Jo = 0.2 and Ko = 0.1 
with increasing I = 10~ 6 , 10" 4 , and 10" 2 . The three dif- 
ferent characteristic energies corresponding to J, K and I 
rapidly merge together as Io increases. 
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